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Abstract 


Total and differential cross sections for high energy and small momentum transfer 
elastic hadron-hadron scattering are studied in QCD using a functional integral 
approach. The hadronic amplitudes are governed by vacuum expectation values of 
lightlike Wegner-Wilson loops, for which a matrix cumulant expansion is derived. 
The cumulants are evaluated within the framework of the Minkowskian version 
of the model of the stochastic vacuum. Using the second cumulant, we calculate 
elastic differential cross sections for hadron-hadron scattering. The agreement with 
experimental data is good. 
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1 Introduction 


In this article we will discuss elastic scattering of hadrons at high centre of mass energy 
a/s (a/s ^ 20 GeV) and low momentum transfer squared f(say |f|^0(l GeV^)). Because 
of the small momentum transfer, such reactions are governed by soft, nonperturbative 
interactions. Experimental data show a rise for the total cross sections of all hadronic 
reactions with increasing centre of mass (c.m.) energy, starting at about a/s = 10 GeV. 
Donnachie and Landshoff (DL) showed that this rise can be described phenomenolog¬ 
ically in terms of Regge theory |Q] by pomeron exchange. The DL pomeron couples like 
a G = -|-1 “photon” to single quarks in the hadrons. The transition from the quark to 
the hadron level leads then to the additive quark rule [Q. Donnachie and Landshoff htted 
the rise of all hadronic cross sections with one small power of s, indicating that there is 
a universal mechanism which governs this kind of reactions. There is also a lot of data 
available for elastic differential cross sections at different c.m. energies, mainly for pp and 
pp scattering [SS0, but also for vrp and Kp scattering [Q . Surely the mechanism which 
governs the elastic amplitude in the forward direction should also control the elastic dif¬ 
ferential cross section (dcr/dt) for sufficiently small |t|. Indeed, pomeron exchange is also 
able to describe da/dt 0. There are many other proposals and methods to describe such 


hadronic reactions, from perturbative held theoretic calculations topological expan¬ 
sions and strings []^, valons 0, the work of Gheng and Wu on high energy behaviour 
in held theories based on perturbation theory [jT^ , to “geometrical” models, which invoke 
global phenomenological properties of hadrons like their “blackness” 

“pomeron physics” we refer to [|I5|. 


For a review of 


A new ehort towards a microscopic description of high energy soft hadronic reactions 
was made in [jl^. In an abelian gluon model the pomeron properties were related to 
nonperturbative properties of the vacuum like the gluon condensate 0 and a “vacuum 
correlation length” a 0]. In [|T^ these ideas were generalised for QCD. It was shown 
there that the amplitude for qq-scattering at high energies is governed by the correla¬ 
tion function of two lightlike Wegner-Wilson lines. Using this formalism a description of 
hadron-hadron scattering was developed in |^, where the hadronic amplitudes are 
calculated from correlation functions of lightlike Wegner-Wilson loops. These correlation 
functions are evaluated in the model of the stochastic vacuum |^, applied in Minkowski 
space after an analytic continuation from Euklidean space. There are by now many other 
applications of this technique, for example in the description of exclusive vector meson 
Related techniques have been used in for dealing with hard diffrac- 


production 

tive processes in deep inelastic lepton-nucleon scattering observed at HERA 


The goal of our paper is, to use and further develop the description of high energy 

In Sect. 2 we collect the formulae for 


diffractive hadron scattering given in Hl^, ^ 


the hadronic scattering amplitudes as derived there. We begin Sect. 3 with a summary of 
the model of the stochastic vacuum (MSV) in its Minkowski version. In the second part 
of Sect. 3 we calculate the correlation function of two Wegner-Wilson loops, the main 
ingredient of the meson-meson scattering amplitude, using a matrix cumulant expansion 
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and the MSV. A question of interest is whether or not the constituents of the baryons 
prefer quark-diquark like conhgurations where two quarks are close to each other on a 
scale given by the proton radius. In |^, ^ strong arguments for the quark-diquark 


picture were given where baryons do act in a hrst approximation as colour dipoles in 
the same way as mesons. In Sect. 4 we present our results for the pp and pp elastic 
differential cross section da/dt treating the baryons as such colour dipoles. In Sect. 5 we 
discuss meson-meson and meson-baryon scattering considering baryons again as colour 
dipoles. Treating meson-baryon and baryon-baryon scattering for general three-quark 
baryon conhgurations along these lines is more complicated and we defer this to another 
publication. Our conclusions are drawn in Sect. 6. 


2 The hadronic amplitudes 


Consider elastic scattering of two hadrons hi, ^2 in the c.m. system at high energies and 
small momentum transfer 


hi (Pi) + h2(P2) ^ hi(P3) + h2(P4). 


( 1 ) 


Now look at this reaction with a microscope. We have to choose an appropriate resolu¬ 
tion in order to extract the essential features of the reaction, but not resolve unimportant 


details. In [|I^ this resolution was estimated, based on the uncertainty relation, with the 
following conclusions. Over a time interval to ~ 2 fm around the “nominal” interaction 
time (i) the parton state of the hadrons does not change qualitatively, i.e. parton annihi¬ 
lation and production processes can be neglected, (ii) partons travel in essence on straight 
lightlike world lines and (iii) the partons undergo “soft” elastic scattering governed by 
non-perturbative gluon dynamics. Using this approach hadronic amplitudes for high en¬ 
ergies were derived in [^, The result for meson-meson scattering where mesons are 
represented as qq wave packets is 

% = 6f, + *(27r)"5(P3 + P 4 - Pi - P2)Tf„ 


Tfi = (-2is) J d%TQ:Kp{iqTbT) J d^xrd'^yTW^iixT) 


w’f(hT,XT)w’:‘(-hT,vr) - 1 ;^ 


( 2 ) 


Here the assumption is made that the q and q share the longitudinal momentum of the 
meson roughly in equal proportions. The interpretation of (^) and the symbols occurring 
there is as follows. The scattering amplitude is obtained by hrst considering the scattering 
of quarks and antiquarks on a hxed gluon potential and then summing over all gluon 
potentials by path integration, indicated with the brackets ( )g- Travelling through a 
gluon potential the quarks and antiquarks pick up non-abelian phase factors. To ensure 
gauge invariance the phase factors for q and q from the same meson are joined at large 
positive and negative times, yielding lightlike Wegner-Wilson loops W± which are dehned 
as 

vr" = brl/(c±), (3) 
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(4) 


V(Ci)^F expl-ig dx^ G‘{x)^]. 

Here V{C±) are non-abelian phase factors (connectors) along cut loops C± as sketched 
in Fig. 1^. The trace in corresponds in the usual way to the closure of the loop. The 



Figure 1: The light-like Wegner-Wilson cut loops in Minkowski space time, C±, consisting 
of two light like lines in the hyperplanes x=p = 0 and connecting pieces at infinity. The 
loops are cut open at one corner, and A+, respectively. In transverse space the centres 
of the loops are at ± 67 -/ 2 , the vectors from the antiguark to the guark lines are Xt and 
yrp respectively. 

transverse separation between the centres of the loops is given by 67 -. The vectors xt 
and yrp give the extensions and orientations of the loops in transverse space. The path 
integration correlates these loops and so causes the interaction. The resulting loop-loop 
correlation function has to be integrated over all extensions and orientations of the loops 
in transverse space with a measure given by the meson’s overlap functions w^i and 
for which one has to make a suitable ansatz. In order to obtain the hadronic amplitude a 
Fourier transform over the impact parameter 67 - has to be done hnally. In the following 
we call (|^) the meson-meson amplitude. As discussed in the introduction, we can use it 
also to describe meson-baryon and baryon-baryon scattering treating baryons as colour 
dipoles in the quark-diquark picture. 


3 Evaluation of the meson-meson scattering ampli¬ 
tude 

In this section we will use the model of the stochastic vacuum (MSV) to perform the 
functional integral in (H) in an approximate way. 
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First we give a short summary of the relevant properties of the MSV. For a detailed 
discussion see 


21 


I . The most important ingredient is a special ansatz for the correla¬ 
tion function of two parallel transported gluon held strength tensors, shifted to a common 
reference point o along the curves and 


Here the right-hand side of @ depends on X\,X 2 and The reference point o can 

be freely shifted on the curve C 12 = Cx^ + Cx 2 where Cx 2 is the oppositely oriented curve 
Cx 2 - The correlation function is proportional to <5“^ due to colour conservation. Now 
the MSV makes the assumption, that F^j^pa is independent of the connecting curve C 12 . 
Then Poincare and parity invariance require the correlator to be of the following form 
{z = xi- X2) 


Fpupcr{z) ^^G2{{g/J,pgua 9p,a9up)\f^ D{^Z ) -|- (1 k)Di{^Z )] 


AD (z^) 

H“ fip ^pa “1“ ^p^pQua ^a^p9i/p)i^ 


Here G 2 is the gluon condensate, D and Di are invariant functions normalized to hl(0) = 
Ili(O) = 1. For spacelike separations they are assumed to fall off rapidly on a length scale 
given by the correlation length a ~ 0.3 fm. The Fourier decomposition of the correlation 
functions is given by 


We follow the authors of 121 


and take as ansatz for D and Di 


b{k^) 

biib) 


27{2np ik^ 

(8a)2 {k'^ — X~‘^ + ieP^ 

2 27{27ip i 

3 (8a)^ (/c^ — A“2 -I- ieP 


( 8 ) 


The constant A is given by A = Sa/Sn and k is a parameter related to the non-abelian 
character of the correlator \21, |^ . 

The Euclidean version of the correlator (^) has been investigated in lattice QCD . 
The ansatz (|^, (H) gives a good description of the nonperturbative part of the correlator 
in comparison to the data from the measurements in quenched QCD and from a £t one 
finds the following ranges for the parameters G 2 ,a, k p3[| : 


nG 2 a^ = 0.39 to 0.41, 

K, = 0.80 to 0.89, 

a = 0.33 to 0.37 fm. (9) 
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From the lattice data one obtains directly the dimensionless quantities k and KG 2 a^- 
Their uncertainty is due to statistical errors of the lattice data and to variations of the 
£t range chosen. The values for a depend also on ALattice which introduces some further 
uncertainty. 

It was shown in [Q that k 7 ^ 0 is crucial for deriving confinement in the framework of 
the MSV. As was found for the total cross sections in and as we will find again here 
a value k 7 ^ 0 is also crucial to describe the experimental data on high energy scattering 
for da/dt in the framework of our model. 


3.1 A cumulant expansion of the colour dipole correlation func¬ 
tion 

Here we calculate the colour dipole correlation function {W^W^)g introduced in (|^). 
The strategy is the following. We transform the line integrals and into a surface 



Figure 2: The curves C+ and C- along which the connectors in W± are taken. The 
mantle of the pyramid with apex at the origin of the coordinate system and boundary 
C+(C_) is P+{P-), the basis surface S+{S-). 


integral using the non-abelian Stokes theorem ||^. Following we choose as surface 


the mantle of the double pyramid P = P+ + P_ which has C+ + C- as boundary and o 
as apex (Fig. H). 

To give the details, let P± be the pyramid surfaces excluding the base surfaces S± and 
cut from o to A±. Then the boundaries of P± 


are 


9P+ = oA+ + (7+ + A+o, 

dP- = oA_+C- + A_o. (10) 


The methods developed for the non-abelian Stokes theorem as explained in [ 2 f] allow us 
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to write the line integrals of @ as 

V{C^) = V{A^o)V{P^)V{oA^), 

V{C_) = V{A_o)V{P_)V{oA_). (11) 

Here V{oA±), V{A±o) are connectors along the straight lines from A± to o and o to A±, 
respectively. They satisfy 


V{oA+)V{A+o) = 11 

V{oA_)V{A_o) = 11 (12) 

The matrices V{P±) in ([TlD are surface-ordered exponentials of held strength tensors G 
parallel transported to o: 


V{P+) = P exp 
1/(P_) = P exp 




A“ 


2 2 


2 jp_ ^ 2 

Inserting (pT]) in (2-3), using (p!2D and the cyclicity of the trace, we get 


(w'f (hr, Xt) Vt)}^ = 

= (l|trr(Cn]l[trV^(C_])^ 

= (5|trr(F+)]i|trV-(P_])^. 


(13) 


(14) 


The main idea is now to interpret the product of the two traces (tr) over 3x3 matrices in 
(P^ as one trace (Tr 2 ) acting in the 9-dimensional tensor product space carrying the prod¬ 
uct of two SU(3) quark representations. Using the dehnition of the matrix multiplication 
in the product space giving e.g. 


(A“0 1)(A^®1) = A“A^®1, 

(A“(8)l)(l® A^) = A“®A^ 

and of path ordering we get immediately 

(H-'f lTV2(Pexp|-| 1 )] 

P«xp[-y ^ <i(T»“'G;„(l® y)| (16) 


The two exponentials in (16) commute because the two matrix structures in the exponents 
do. Introducing a total shifted held strength tensor Gt as 


GGx) 


GjAo,pCx){^®G for xeP+ 

for xeP- 


(17) 
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we can rewrite the two exponentials in ([T6|) as one exponential defined in the direct 
prodnct space. In this way we get from (|I^) a path-ordered integral over the donble 
pyramid mantle P = P+ -|- P_: 


W^W 


m\ 


- / 


^ = ^Tr 2 (P exp[-i| da{x) Gt{x)] 


G 


(18) 


Here and in the following we snpress the Lorentz indices if there is no confnsion. Note that 
the path orderings on P+ and P_ do not interfer with each other. Thus the path-ordering 
on P can for instance be chosen such that all points of P+ are “later” than all points of 
P_. 


For the expectation value of the single surface ordered exponential (jlSD we can make a 
cumulant expansion [^, ^ . In our case we use a matrix cumulant expansion as explained 
in (2.41) of |2^ (cf. also [061): 


9 f '' ^ ^ 9 f 

(Pexp[-i^ / dcr(a:)Gt(a;)]) =exp[ -r/ da{xi) ■ ■ ■ da{xn) Kn{xi, ..,Xn)]. 

' z Jp ' G 77,1 Z J 

(19) 

Here the cumulants Kn are functional integrals over products of the non-commuting ma¬ 
trices Gt of (O)- Thus one has to be careful with their ordering. The cumulants up to 
n = 2 are 


Ki{x) = {Gt{o,x;G^))(., 

K2{Xi,X 2) = {P[Gt{o,Xi;G^^)Gt{o,X2;Ga:2)])Q- 

Gt{o,xi]Gxi))Q{Gt{o,X2',Gx2))G + (1 2) V 


( 20 ) 


Note that the Gt have Lorentz indices and are matrices in colour space as shown in 
(0- The functional integral indicated by ( )g in (|^) involves only the field strength 
components thus also the cumulants Ki, K 2 , ... still carry Lorentz and colour indices. 
The fact that there is no colour direction preferred in the vacuum requires Ki to vanish. 
Neglecting cumulants higher than n = 2 we get for ( 0 ) 

= ixtj exp(C2(a:T, Vt, 6t)), 

2 

G2ixT,yT,bT) = J^da{xi) J^da{x2) (P{Gt{o,xi;Gx^)Gt{o,X2;Gx2)))^-i‘21) 
where (^2 is a 9 x 9 matrix, invariant under SU(3) colour rotations. 


3.2 Calculation of the second cumulant term using the MSV 

Now we use the MSV with the ansatz (0),(^ for the invariant functions D and Pi to 


calculate G 2 of (0ID , which can be split into three contributions: 


py pyP+P+ _i_ G^P-P- I G^P+P- 

C 2 — 0/2 -r G’2 <^2 


(22) 







p p p p 

In each of C' 2 ’^ and C 2 ~ ~ both xi and X 2 move on the same snrface, P+, P_, re- 

p.p 

spectively, and so we have to pay attention to the snrface ordering. In € 2 ^ ~ one point 
moves on P+, the other one on P_. In this case it follows from (p!7D that the shifted held 

p p 

strengths commute and the surface ordering is irrelevant for € 2 ^ ~■ 

Now we show that vanishes. To see this we follow the argumentation of 

and transform the two surface integrals over the pyramid mantle P+ in (|^) into surface 
integrals over 5*+ and integrals over the volume Id+ enclosed by P+ and 5*+. The integrals 
over 5+ vanish due to the Lorentz structure of the surface elements together with 
the ansatz (^. The integrals over I/+ are roughly speaking sums of integrals over surfaces 

p p 

parallel to 5"+ and so they vanish for the same reason. In a similar way we see that 6*2 “ “ 
vanishes. 

To calculate € 2 ^ “we use the same method. In this case neither the integrals over S± 
nor the ones over V± vanish. The integrations in light-like directions can be done analyt¬ 
ically. Using then also the ordinary Gauss theorem we hnd with (^-(|^ that everything 
reduces to line integrals over the vectors Vyi {i = q, q) running from the apex o to 
the position of the quarks and antiquarks in transverse space (Fig. 


A^ 

C 2 {xt, Vt, bn) = y Y^“*) 


x{xt, Vt, br) = 


24 “ 1 “ f'yq) 

I{Xxq-i f^yq) ~ ^i^xq, ^'yq) } j 


I(r„ r,) = * [' dv, f dv, r r, - v, r,) 

Jo Jo J -00 (2 vr)^ 


{nVyr^Di-kl) + {1 - k) {kTry){kTr^) P((-fc|)}, 




(23) 


As we can see we hnally need the correlator functions D,Di of (^) for space-like momenta 
only. This means that the result (^) involves the correlation functions D{z^) and Di{z^) 
only for spacelike z, where they are as in Euclidean space time. With (H) for the functions 
P, Pi we get 


K^x^ry) = {Vyr^) ^ dn((- 


PjLZlAfK 


\-y---AfK2{ 


A 

\r — vr. 


) + 


A 


A 


')) 




(24) 


where 7 ^ 2,3 are the modihed Bessel functions. Note that y is a real function. 
Inserting (E^ in (|^) we obtain 
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Figure 3: The projection of the pyramid surfaces P± + S± into transverse space. The 
vectors r^^q^q) and ry^q(^q) point from the origin o along the transverse projections of the 
pyramid mantles P± to the projections of the guark and antiguark lines of C±. xt and 
yrp point from the antiguark to the guark lines. 


Xt) Vt))^ = 

- Trs exp[(y 0 —) {-i) x{xt, Vt, br) 


To evaluate the trace in (p5[) we introduce two projectors Ps and Pa 


(-Ps)(ai02)(/3i/32) 2 *^“2/32 T daig2^oi2l3i) J 

(■Pa)(oia2)(/3i/32) 2 *^“2/32 <^«i/32*^02/3i) 


(25) 


(26) 


which act in the direct product space of two SU(3) quark representations projecting onto 
the subspaces carrying the irreducible representations. The decomposition is: 3®3 = 6©3. 
Using the identity: 


y 


0 


A“ 

y 




(27) 


together with the projector properties of Pg and Pa and Tr 2 Pg = Q and Tr 2 Pa = 3 we can 
immediately calculate the trace in ([25|) and so the colour dipole correlation function: 


A^ 

Tr2exp [(y 0 Y)(-f x)] 

= Tr2exp[(^P, - ^Pa)(-fx)] 

= Tr2[P,e-*5^ + P,pi^] 

= 6e"*i^ + 3e*i^. (28) 
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3.3 The meson-meson amplitude 


Inserting (|^) in (y) and using the Gaussian shaped mesonic overlap functions from pi 
our hnal result for the meson-meson scattering amplitude reads 


Tfi = (2is) (2 7r) / 

( 

JM,M{b) = — (f'XT 


dhhJQ{^\h)JM,M{b), 
1 


2 TT 


exp 




(fyT 


^H2 


exp 


2 1 12 

2 cos(-x(a;r, Vt^ &r)) + ^ cos(-x(a;T, 2/r, ^t) ) + 

2 1 12 
- 2sin(-x(a;T,yT,^T)) + - sin(-x(a::T, 2/r, &t) ) 


Vt 




(29) 


Here Jq is the zeroth-order Bessel function. Due to ( p8D we get from the original matrix 
valued exponential (^) a sum of two c-number valued exponentials which we have written 
in terms of trigonometric functions. Different hadrons are distinguished in (p9|) only 
through their strong interaction extension parameters Sh which should be of the order of 
the electromagnetic hadron radii. 

Assuming |x| -C 1 and expanding in (P^) to the order O(x^) gives the result 

of pg. 


Tf, = {2zs) ( 27 r) r dbbM^\b)f^]^{b), 


— j d'^XT 


2^rS^ 


exp 


Xn 




d'^VT 




exp 


yl 


2SJ. 
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(x[xT,yTM)?, 


(|xl«l). (30) 


But the integral JM,M{b) in (^) is dominated by a region in xt^Vt where |x| <C 1 only 
for larger values of b, say 6 > 4a (see section 4). So using (|30[) JuM^b) can be calculated 
reliably only for larger values of b. Nevertheless with (|5I1[ ) the total cross section and the 
slope parameter at t = 0, where one needs only the hrst and third moments of JM,Mib), 
were calculated in in a satisfactory way. 

Coming back to our expression (^) for the elastic scattering amplitude we show next 
that this amplitude is purely imaginary. Since y is real, any real part of (^) would have 
to come from the sine-terms. Now the overlap functions in (P9[ ) are invariant under the 
replacements: xt —>■ —xt and yrp —*• —yrp, respectively, but 


x(-a^T, 1/r, br) = x{xt, -!/t, ^t) = -x{xt, yr, br)- (31) 

This is easily seen from (^) since xt —»■ —xt means the replacement r^q Vxq and 
yrp —>• —yrp the replacement Vyq ^ Vyq. Thus integrating over xt and y^r averages out 
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the real part in and we get 

Tfi = {2is) {2 7 i) f dbbJo{\/\t\b)JM,Mib), 


JM,M{b) — — d?XT 








d^UT 


2¥^"^P 


yl 


2^l.„ 


2 1 12 
g cos{-x{xt, Vt, &t) ) + g cos(-x(a::r, Vt, &r) ) - 1 


( 32 ) 


As a consequence our meson-meson amplitude is invariant under the replacement of one 
hadron by its antihadron. The exchange of all partons by their antipartons for a given 
parton conhguration turns around the loop direction. This results in a change of sign 
of X and so does not affect the amplitude (^). In our approximations, we get only 
charge conjugation C = -|-1 (pomeron) exchange and no C = —1 (odderon) exchange 
contributions to the amplitude. A real part of the amplitude and C = —1 exchange 
contributions could arise from higher cumulants in (pT]). 

We discuss now the constraints on the elastic amplitude implied by the partial wave 
unitarity (see for instance |^). The partial wave expansion for the T matrix element for 
spin-zero mesons reads 


T{s,t) = ^{21 + l)Pi{cosy)ai{s), 

cm ;=o 


(33) 


0<vi<^- (34) 


Here Pi are the Legendre polynomials, Pcm is the cm momentum, ^ the cm scattering 
angle, and 6i,rii are the phase shifts and inelasticities, respectively. 

At high energies we get from (|33|) with b = {21 + 1)/s/s\ 

T{x,t) = 87is J dbbJQ{b^J\t\)al{^/s). 

Comparison with (^) gives 

JMM{b) = + 1 

Thus, the partial wave unitarity requires 

\JMM{b) — 1| < 1 (37) 


2l+l=b^ 


(35) 


(36) 


and this is always satished for our amplitude (|3^) since x is real and the prohle functions 
are normalised to one. 
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4 Proton-proton and proton-antiproton scattering in 
the quark-diquark picture 


A lot of experimental data is available for elastic differential cross sections of proton- 
proton (pp) scattering up to ^/s = 63 GeV and proton-antiproton (pp) scattering up to 
Tevatron energies of y/s = 1800 GeV ii 0- In this section we will compare the high 
energy data in the range y/s > 23 GeV to the results of our calculations making the 
assumption that the proton has a quark-diquark structure. Thus we can use the formulae 
for the meson-meson amplitude as presented in Sect. 3. 

Our starting point is (^) depending on 4 parameters: the QGD vacuum parameters 
G 2 , n and a and the proton extension parameter Shi = Sh 2 = Sp. The vacuum parame¬ 
ters are surely energy and process independent, the extension parameter Sp will be allowed 
to vary with energy. The numerical calculations using (|3^ are too lengthy to attempt a 
“best £t” of these parameters from the data. We adopted the following procedure instead. 
In the SVM with the ansatz (0|D for the functions D, Di the string tension p is given by 


P 


TT^kG2 
36 Jo 
327iKG2a‘^ 


dZ^D{-Z^) 


81 


(38) 


Typical values for p extracted from phenomenology (cf. e.g. [0) are p = (420±20 MeV) 


In the following we will express G 2 through p and a using (|38|) . From previous work 
and the lattice measurements discussed in Sect. 3 we expect for the correlation length 
0.30 fm ^ a ^ 0.37 fm, for n ^ 0.75 and for the proton extension parameter at y/s = 23 
GeV Sp ~ 0.86 corresponding to the electromagnetic proton radius. Now we considered 
again a/s = 23 GeV and started our numerical investigations using (|3^) , varying the 
parameters p, a, k, Sp aroung the values indicated above. We calculated 


— =-iTfiP 

dt 167r s^ 

and the total cross section, using the optical theorem which reads for s 3> m^\ 

aripp) = ^Im(T/i) 


(39) 


(40) 


and compared to experiment. The experimental data on the total pp and pp cross sections 
is very well described by the DL fit [0, |l|. We are only interested in the pomeron exchange 
part here. Therefore we take as “experimental” input the pomeron part of (Tt{pp) in the 
DL parametrization p[: 


aripp) 


exp 



0.0808 

mb. 


(41) 


We imposed as constraint that our amplitude reproduced ( 0 ) exactly at yfs = 23 GeV. 
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With this procedure we found quite a satisfactory description 
shown in Fig. 5 for the following values of our parameters: 

of the data for da/dt as 

p = (435 MeV)^ 

(42) 

a = 0.32 fm, 

(43) 

K = 0.74, 

(44) 

Sp{s = (23 GeV)^) = 0.87 fm. 

(45) 

From (H), (H), (H) we get 


G 2 = (529 MeV)^ 

KG2a^ = 0.40. 

(46) 


Varying the parameters away from the values did not lead to improvements. Also, 

the values (|4^ - |45|) are well within the range obtained in pl|, |30| and quite compatible with 
the lattice results (^. In the following we will thus £x the vacuum parameters p, a, k, to 
their values 

Now we consider aT and da/dt at higher energies ^/s^ where we have then only one 
free parameter Sp{s) left. We fix Sp{s) again by requiring that our model reproduces the 
experimental value for the pomeron part of cxipp) according to (^I[) . 

p7| and £t our calculated values of (Jt{.PP) from 
We get a good description in the range 
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On the other hand we follow 
(^2]), (l40|) as shown in Fig. | to a power of Sp 



Figure 4: The total cross section ax versus Sp. The solid points are the calculated values, 
the solid curve corresponds to the power fit. 


2.5 a < Sp < 4.0 a with 


MSp) 


cal 


50.384 {Spf^^^ mb. 


(47) 
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Equating (|40|) to (|4^) we get Sp as function of s: 


Sp{s) 


0.700 



0.034 

fm. 


(48) 


This leads to Sp = 0.86, 0.87, 0.93, 1.07 and 1.17 fm for ^/s = 20, 23, 63, 546 and 
1800 GeV, respectively. For comparison, the mean squared charge radius of the proton as 






Figure 5: Differential elastic cross sections for c.m. energies ^/s = 23,63,546 and 1800 
GeV. This corresponds to proton extension parameters Sp = 0.87, 0.93, 1.07 and 1.17 fm. 
The data at a/s = 23 and 63 GeV are from the ISR experiment p . The pp scattering 
data at a/s = 546 GeV are from the and the data at a/s = 1800 GeV from 0]. 


determined from Lamb-shift measurements 


is rp = 0.89±0.014fm. Now everything is 

we 


fixed and we can compute the elastic scattering amplitude and ap from (p2D. In Fig. 
present our results. A first observation is that for all energies the calculated differential 
distributions follow the experimental data quite well over many orders of magnitude. The 
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fact that this is true up to y/s = 1800 GeV supports the description of the s-dependence 
by a s-dependent extension parameter Sp{s). 

For small |t| ~ 0.25 GeV^ we get a change of slope (Fig. |^). Splitting the integrand 



Itl / GeV^ 


Figure 6: Differential elastic cross sections calculated for the values of G 2 = (529MeV)'^, 
a = 0.32 fm and Sp = 0.87 fm using as integrand of the amplitude (|3^) |[cos( | y) — 1] 
(solid line) and |[cos( | y) — 1] (dashed line) respectively. 


of (|3^) into two contributions 


o cos(-x) 


o cos(-x) 


1 = 


r 2 , 


2 


1 X 1 

cos( - X) 

-1 

^ 3 

cos( 

3 2 - 


(49) 


we hnd that for |t| ^ 0.25 GeV the hrst term dominates and for |t| ^ 0.25 GeV the second 
one dominates. Such a change of slope is indeed reported by experiments . 


For all energies the imaginary part of our amplitude changes sign at some t < 0. Due 
to the absence of a real part in (^) the calculated differential cross sections have a zero 
there. This causes an inhnitely deep dip in our t-distributions. We expect this dip to be 
at least partly hlled up once we change to more general quark conhgurations and include 
higher cumulant terms. The point at which the zero occurs in our calculation moves to 
smaller values of |t| with increasing energy and is always in the region where experiments 
see a marked structure: At lower energies there is a dip in pp, and a shoulder in pp 
scattering, respectively. At the highest energies only pp data is available and one hnds 
a shoulder. Thus our model produces structure in the t-distributions at the right place. 
But we should insert the warning that these dips occur at |t| ~ 1 GeV^ where one would 
expect also perturbative effects to play a signihcant role. (cf. P0|, [p[| ). 

Of course, our model does not give a perfect £t to the data. At all energies our 
calculated curves are somewhat too steep at very small |t|. Our amplitude is purely 
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imaginary and thus does not satisfy the relation between the phase and the s-dependence 
required by analyticity and Regge theory [^. Also our da/dt is the same for pp and 
pp scattering, whereas experimentally these differ markedly in the dip region. This was 
nicely explained theoretically in as an interference of single and double pomeron and 
three-gluon exchange. We will have to see if higher cumulant terms and/or a departure 
from the strict quark-diquark picture of the proton will lead us to an improvement on 
these points in our model. 

In Fig. 1^ we show our prediction for da/dt in pp scattering at \/s = 14 TeV corre¬ 
sponding to the LHC energy. The total cross section is again assumed to be given by the 



Figure 7: Prediction of the differential elastic cross section for the LHC energie ^/s = 
14 TeV. 


DL parameterisation (^). This was used as constraint to fix the extension parameter Sp 
according to Fig. 4 with the result 


= (14 TeV)^) = 1.34 fm. 


(50) 


But the shape of da/dt is a prediction of our model and we are looking forward to the 
corresponding experimental data. 

We will show next that our results for da/dt depend crucially on the string tension 
p 7 ^ 0, i.e. on the confinement features of QCD. For this we plot in Fig. ^ da/dt for 
a/s = 23 GeV calculated for the same values of G 2 = (529MeV)^, a = 0.32 fm and 
Sp = 0.87 fm, but for different values of k. For the purely non-abelian case k = 1 the fall 
of da/dt is too steep, for the purely abelian, the non-confining case k = 0, the fall is much 
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Figure 8: The variation of the differential elastic cross section calculated for the same 
values of G 2 = (529MeV)^, a = 0.32 fm and Sp = 0.87 fm; (a) for k = 0.74 (solid line), 
n = 1 (dashed line) and k = 0 (dash-dotted line); (b) for n = 0.74 (solid line), n = 0.90 
(dash-dotted line), n = 0.80 (dashed line) and n = 0.60 (dotted line) together with the 
experimental data (solid points) at ^/s = 23 GeV. From (b) one can see, that a value of 
K around 0.74 is favoured by the data. 

too slow with increasing \t\. The correct fall and dip position is obtained for k = 0.74, 
not far from the value determined from the lattice calculations (H). 

To explain this we plot the profile function JM,M{b) multiplied with 6 for k = 0 and 
K = 1 versus b in Fig. p. As we can see, the two curves differ noticeably. For k = 1 we 
get a single maxima whereas for k = 0 we get two maxima. This can be understood by 
looking at the dependence of y on 6 for hxed transverse vectors Xt and in the two 
cases. As we see from (p3D y is a sum of four terms (replace here q by the diquark qq) 
corresponding to q — q, q — qq etc. interactions. For k = 0 the function /(r^,, Vy) in 
depends only on the difference ry — r^ of the parton positions, but for k = 1 lir^, Vy) gets 
contributions from all “strings” spanned between o and the g’s and gg’s positions ( Fig. 
|T^). We get large contributions to I if the arguments of the Bessel functions K 2^3 in (^) 
are small. For k = 0 this means that the partons (g, qq) have to be close in transverse 
space, but for k = 1 it is only required that one parton be close to the string of the other 
parton. To see the consequences of this, we take as an example \xt\ = |?/t| = 4 and 
Xt, Vt, bx nearly parallel to each other (Fig. m- For K = 1 we have the following 
situation. For small and medium values of b the terms I corresponding to the q^ — qy 
and qqx — qqy interaction dominate in the sum for y in (^). These functions I decrease 
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Figure 9: The dimensionless profile function JuMip) times b versus b for G 2 = 

(529MeV)"^, a = 0.32 fm and Sp = 0.87 fm in the cases k = 1 (solid line) and k = 0 
(dashed line). 

with increasing b nearly monotonously corresponding to moving away from the string 
o — Qy and QQy away from the string of o — qq^. But the interaction of qy with the string 
o — qx and qq^ with the string o — qqy remain dominant. Compared to them the terms I 
corresponding to the qq^ — qy which enters with negative sign in ( 1 ^) stays smaller and 
also the function y decreases nearly monotonously with increasing b. In summary: for 
K = 1 the function y is sizable and practically always of the same sign as long as partons 
are close to strings of other partons, i.e. as long as projection of the two loops overlap. 
This picture, smeared out due to the integrations over the transverse vectors, is reflected 
in the shape of {bJM,M{b)) for k = 1. 

For K = 0 the function y has a maximum at 6 = 0, where the quarks and the diquarks 
of the two protons are closest to each other in transverse space. Then y decreases rapidly 
for increasing 6 on a scale given by the correlation length a, passes zero at some value 60 
and approaches a minimum. The latter corresponds to the situation where the quark of 
one proton is very close to the diquark of the other proton. Since there are no strings 
in this case all interactions except the one between qq^ and qy are then negligeable and 
the latter enters with opposite sign to the qx,qy and qqx,qqy interaction in (^3]). For 
b increasing further the transverse distance between all partons of the dipoles increase 
and y goes to zero. At bo the expression in square brackets in (^) vanishes. Thus we 
expect - after smearing out through the integrations - to see a minimum in b-TMuib) 
around bo and this is indeed seen in Fig. p. Thus, through this sign change of y we 
understand the two maxima in the shape of 6 for k = 0 , and also that for k = 0 

the hrst maximum in bJMM{b) occurs for a smaller value of b than for k = 1. After the 
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Figure 10: (a) The transverse positions of the quarks q^, qy and diquarks qq^, qqy for 
parallel loops. The dashed lines represent the “strings” spanned between o and the q’s 
and qq’s positions, (b) The eorrelation function x{b) versus b for \xt\ = ly^l = 4a and 
Xt, Ut) br parallel to each other in the cases k = 1 (solid line) and k, = 0 (dashed line). 
The curves correspond to loop configurations sketched in (a), namely to parallel loops lying 
on each other for 6 = 0 and moving for increasing h in +xt direction and —yx direction 
respectively. 




Fourier-Bessel transformation in P^ ) this translates immediately in the slower decrease 
of da/dt for k = 0 compared to k = 1. The dip in da/dt is generated by a cancellation 
of positive and negative contributions from bJMM{b)Jo{ba^J]t\). For k = 1 and k = 0.74 


the function hTMuib) is “smooth” and the oscillating Bessel function brings the integral 
(^2]) “easily” to zero for some |f|. For n = 0 the “oscillation” of h-TMuip) together with 
the oscillation of can produce a dip only at much higher |t|. From this point of 

view the structure of the pp scattering amplitude gives us direct information that a string 
formation between quarks in QCD is essential also here. 

We can now also understand how the shrinkage of our calculated differential cross 
sections (see Fig. |^) for increasing c.m. energies is caused: We found above for k = 0.74 
that the shape of da/dt is controlled by string scattering. We assumed Sp to increase 
with s which means increasing lengths of the string sizes and thus still overlap for larger 
values of the impact parameter 6. This in turn translates into steeper forward peaks of 
the scattering amplitude and so into steeper elastic differential cross sections at higher 
energies. 
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Figure 11: 

For ^/s = 23 GeV (a) the profile functions JMM{b) (solid line) and (dashed line) 

as a function ofb. (b) The profile functions JMM{b) (solid line) and J^Muib) (dashed line) 
multiplied with b versus b. 


In Fig. I^we compare our profile function JMM{b) from (^) with the corresponding 
one, J^MAdib) of (10)) where only the term of order is kept. We see from Fig. 0 that 
this can only be a good approximation for larger values of b and indeed, the two functions 
Jmm and Jmm are quite similar for b ^ 2a, but are rather different for b ^ 2a. But for 
the total cross section and the slope parameter at t = 0 which are obtained from integrals 
over bJuMip) and fi^JuMip) this does not make much difference. On the other hand at 
larger and larger |t| the integrals with the Bessel function in (R^, (pOD probes smaller and 
smaller vales of b in JMM{b) and then the differences between Jmm and Jmm show up 
clearly. 

Finally we come back to the partial wave expansion for the T matrix element at 
high energies (^5]) and discuss the case of total absorption which means rji = for the 
inelasticities rji in (PB]). Thus total absorption requires \JMM{b)\ < 1 which is again 
satisfied for our amplitude (^21). At ^/s = 23 GeV this can be read off from Fig. IHa. 


5 Meson-meson and meson-baryon scattering 

5.1 Meson-meson scattering 

Meson-meson scattering is the reaction best suited to apply our formula (^). Unfortu¬ 
nately there is no data available for small t elastic meson-meson-scattering at sufficiently 
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Figure 12: 


Vector meson production in two photon processes at e'^e colliders. 


large c.m. energies. But the reaction 

7(?i) + 7 (^ 2 ) ^ Vi{p^) + 172 (^ 4 ), VI ,2 = (51) 


is nearly as good to extract mesonic differential cross sections as a purely mesonic one and 
it can be studied for instance at LEP [E] in photon-photon processes (Fig. ^). Certainly 
at least part of the reaction (|^) is due to the vector dominance model mechanism: 
The photons fluctuate before the collision into vector mesons which then interact. The 
hadronic interaction can be calculated in the way described in Sect. 3. Thus in essence we 
can regard the reaction (^l|) as elastic scattering of vector mesons. This approximation 
should be best for quasi real photons in (|^ , i.e. for very small virtualities Ig^l and |g||. 

In this spirit we discuss now pp elastic scattering as an example. In our amplitude 


(p2D we have already fixed the vacuum parameters p, k and a in Sect. 4, but we still have 
to choose the p extension parameter Sp and we assume it to be equal to the tt extension 
parameter as determined in the next section 

Indeed, data on p-N total cross sections from photoproduction on nuclei at low energies 
H show (Jt{pN) k, axiT^N). Thus we assume: 


Sp = 0.60 fm 


for \/s = 20 GeV. 


(52) 


Now everything is fixed and we can calculate da/dt for p-p scattering. 

In Fig. |l^ we plot da/dt normalised to 1 at f = 0. Assuming the vector dominance 
model and neglecting contributions from non-diagonal scattering uu —> pp etc., Fig. 
also gives the shape of da/dt for 77 —pp-scattering. It would be very interesting to have 
data for instance from LEP2 to compare with our f-distribution. 


5.2 Differential cross sections for meson-proton scattering 

Here we calculate elastic differential cross sections for 7 r=*"p and K^p scattering at ^/s = 
19.5 GeI7. This is the largest energy for which data for these reactions exist. 
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Figure 13: Prediction for the differential elastic cross section of pp-scattering at y/s = 20 
GeV using Sp = 0.60 fm. 


In our scattering amplitude (^) we have already fixed vacuum parameters p, k, a 


^). We consider the proton as a quark-diquark system with extension parameter Sp and 
use Sp = 0.86 fm from (|48|) . To fix the meson extension parameters we again normalise 
our total cross sections to the pomeron parts in the DL parametrisations of the and 
K^p total cross sections [0. 


a{7r^p)\pom. = 22.0 mb, 

a{K^p)\pom. = (53) 


This leads to 


St, = 0.60 fm, 

Sk = 0.55 fm. (54) 

For comparison the electromagnetic radii are [^: = 0.66±0.01 fm, rx = 0.58±0.04fm. 

Fig. 1^ shows our results for the tt^p and K^p elastic differential cross sections. As we can 

see the experimental data are again reproduced quite well. As in pp scattering the slopes 
of the calculated cross sections are somewhat too large for very small |f|. Furthermore 
all our t-distributions are slightly below the experimental data up to |t| ~ 0.5 GeV. This 
could again be due to with the missing real part in the amplitude (^), a problem which 
might be cured once we go beyond the second cumulant approximation. 

We also note that the c.m. energy y/s = 19.5 GeV considered here is not very high 
and - in Regge language - effects of nonleading trajectories are still sizable and more so 
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Figure 14: Differential elastic cross sections for n^p and K^p scattering at ^/s = 19.5 
GeV. This corresponds to Sp = 0.86 fm, St, = 0.60 fm and Sk = 0.55 fm. The data are 
from p. 


for 'K~p, K~p than for vr+p, K^p scattering. Our calculation does not contain any non¬ 
leading Regge-exchanges and thus should agree better with vr+p, K^p than 7r“p, K~p 
scattering. This is not incompatible with the results shown in Fig. |T^. 

Due to the smaller extension parameters of the mesons our calculated t-distributions 
for meson-proton scattering are flatter than those for proton-proton scattering. Also our 
calculation gives dips only around |t| ~ 2GeV^ and we would not believe our model to be 
reliable at such high |t|-values. 

The squares of the ratios of our extension parameters of mesons and the proton are 


=0.49, (^1 =0.41, (^) =0.84. 


'5 




S', 


S 


K 


s. 


(55) 


The corresponding ratios for the mean squared electromagnetic radii are [^, ^ 


% = 0.55, ^ = 0.42, ^ = 0.77. (56) 

P p 7T 
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The ratios (|55|) follow the trends of (^) but certainly are not equal to them. It is 
particularly noteworthy that in our model the facts that (i) the K^p total cross sections 
are smaller than the ones and (ii) the t distributions for K^p flatter than for vr^p are 
both reproduced quantitatively by a smaller extension parameter Sk compared to In 
the additive quark model on the other hand one has to assume a different cross section 
for the scattering of u, d on u,d and u,d on s quarks which is hard to understand since 
the gluon interaction is flavour-blind. 


6 Conclusion 


In this article we have presented calculations of amplitudes for elastic proton-proton, 
meson-meson, and meson-proton scattering at high energies and small momentum trans¬ 
fer. Our model is based on functional integral techniques and an appropriate eval¬ 
uation of such integrals in the framework of the stochastic vacuum model | 
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comparison with previous work |^, |^, we have now made a cumulant expansion for 
the correlation function of two Wegner-Wilson loops instead of an expansion in terms of 
the number of held strength correlators. We found that this latter expansion can only 
be justihed for medium and large impact parameters, but gives, nevertheless, reasonable 
results for the total cross section and the slope parameter at f = 0. 

As parameters in our model we have the QCD vacuum parameters: the gluon conden¬ 
sate G 2 , the non-abelian parameter k, and the correlation length a. In addition we have 
the s-dependent hadron extension parameters Sh{s). 

We hxed the Sh{s) by requiring the total cross sections to agree with the experimental 
values. In this way we found quite a good description of da/dt for pp,pp, vr^p and K^p 
scattering for vacuum parameters G 2 , n, a or equivalently the string tension p, k and a as 
given in These values are close to the corresponding values obtained from lattice 

data 1^, ^], from previous investigations of high energy scattering ^ and from 
various other low energy phenomena (see for a review). The parameters Sh{s) came 
out close to the known electromagnetic radii of the hadrons for energies ^/s ~ 20 GeV 
and were required to increase slowly with s according to (|48|). 


As in previous work ||^ the fact that axi^Kp) is smaller than axiTip) is related in our 
model to a smaller extension parameter Sk compared to 5^. 

The dip structure in pp scattering around \t\ ^ 1.5 GeV^ was seen to depend crucially 
on the non-abelian character of the gluon held strength correlator, i.e. on k ^ 0. This 
leads to the string formation and the area law for the Wegner-Wilson loop with static 
quarks. In our model these strings play again a crucial role in high energy scattering in 
producing the correct t-distributions. 

At very small values of |t| our model gives a change of slope related to the symmetric 
and antisymmetric combinations of the multi-gluon exchange as explained after (^). 

As stated above one of the main ingredients of our model is a cumulant expansion for 
the correlation function of 2 Wegner-Wilson loops which we truncated after the second 
cumulant. But we can easily see that the main features of our model will remain unchanged 
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if the matrix cumulant expansion (|T^ converges. Then the complete snm of cnmnlants 
in (11^) is a 9 X 9 matrix, invariant nnder SU{3) rotations. We will again obtain a 
decomposition for it as for C 2 in (^) ff in terms of the projectors Pg and Pa with invariant 
fnnctions Xs,a mnltiplying them etc. Of conrse the detailed shape of these fnnctions will 
be different. 

Going back to onr amplitnde (|3^ , we can expand the cosines in powers of x- Maybe 
the x^-term - as kept in the work - conld be interpreted as exchange of two non- 


pertnrbative glnons, as single bare pomeron exchange, the x^-ferm as exchange of fonr 
nonpertnrbative glnons, as donble pomeron exchange, etc. At the moment, however, snch 
an identideation is pnrely speculative. 

In all this work we have treated baryons as quark-diquark systems, i.e. we assumed 2 
quarks of the 3 valence quarks of a baryon to be close together in transverse directions. 
This picture is supported e.g. by the investigation of where it is shown to give 
an explanation for the apparent absence of odderon couplings of the proton at small 
|t|. In future work we plan to investigate meson-baryon and baryon-baryon scattering 
treating baryons as 3-quark systems with arbitrary distances between the quarks. Also 
the dependence of the results on the surfaces spanned into the loops (see Fig. ^ in order 
to apply the non-abelian Stokes theorem will be investigated. 

To summarise: We have presented a model where the t-distributions of elastic hadron- 
hadron scattering are related quantitatively to the parameters of the vacuum in non- 
perturbative QCD. The total cross sections which rise with energy were used to fix the 
energy-dependent effective strong interaction extension parameters of the hadrons. A 
theoretical calculation of these extension parameters remains a challenge. 
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